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Abstract 

We consider the 3-D full Navier-Stokes equations whose the viscosity coefficients and the 
thermal conductivity coefficient depend on the density and the temperature. We prove the 
local existence and uniqueness of the strong solution in a domain fl C K 3 . The initial density 
may vanish in an open set and f2 could be a bounded or unbounded domain. We also prove a 
blow-up criterion for the solution. Finally, we show the blow-up of the smooth solution to the 
compressible Navier-Stokes equations in M n (n > 1) when the initial density has compactly 
support and the initial total momentum is nonzero. 
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1 Introduction 

The motion of a viscous, heat-conducting fluid in a domain Q C M 3 can be described by the 
system of equations, known as the Navier-Stokes equations: 

p t + dw(pu) = 0, (1.1) 

(pu) t + div(pu (g>u)- div(/u(Vu + Vu T )) - V(Adivu) + Vp = pf (1.2) 

(pe) t + div(peu) - div(fi;V6>) + pdivu = ^| Vu + Vu T | 2 + A(divu) 2 + ph, (1.3) 
in (0, +oo) X Q, C (0, +oo) x M 3 , and the initial and boundary conditions: 

(p,0,u)\ t =o = (po, ,u ) in O, (1.4) 

(0, «) = (6 b) 0) on (0, oo) x dn, (1.5) 

(p,9,u)(t,x) -» (p 00 ,^ 00 ^) as |z| -> oo. (1.6) 

In this paper, Q is either a bounded domain in M 3 with smooth boundary or an usual unbounded 
domain such as the whole space M 3 , the half space M 3 ^ (#& = 9°°) and an exterior domain in R 3 
with smooth boundary. The known fields h and / denote a heat source and an external force per 
unit mass, p, 9, e, p and u denote the unknown density, temperature, internal energy, pressure 
and velocity fields of the fluid, respectively. The pressure p and internal energy e are related to 
the density and temperature via the equations of state: 

p = p{p,9), e = e(p,0). 
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In accordance with the basic principles of classical thermodynamics, p and e are interrelated 
through the following relationship: 

d p e = — (p- Od g p). 

In view of classical thermodynamics (see the book |15j). the viscosity coefficient p = p(p, 9), the 
second viscosity coefficient A = X(p,9) and the thermal conductivity coefficient k = n(p,9) are 
required to satisfy the natural restriction 

p > 0, 3A + 2p > and k > 0. 

In this paper, we assume that (e,p, k, p, A) satisfy 



and 



dge(p, 6) > c > 0, < p(p, 9) < P ho(p)(l + 2 ), 


(1.7) 


\d e p(p,e)\< P hi(p)(l + 0), \d pP (p,9)\ < P 2(p)(l + e 2 ), 


(1.8) 


p(p, 0) > p° > 0, 3X(p,9) +2p(p,0) > 0, 0) > k° > 0, 


(1.9) 


p-*do(d e e,d p e,p,\,K)(p,e) GL£.([0,oo) x [0,oo)), 


(1.10) 


d e e,d p e,p,p,X,n £ C^^oo) x [0,oo)), p ,Pi, P 2 G C([0,oo)), 


(1.11) 


9 tf ,9 () p€C a ([0, TO ) ) C([0 ) cx)))) ) 


(1.12) 



for all p,0 > and a constant a € (|, 1). It is easy to see that the ideal flow ( p = Rp6, e = 
C v 9, and C v , k, p, A = constants) satisfies the above assumptions with po = pi = pz and P 2 = 1. 
The above assumptions are motivated by the facts in [21 [23] where it is pointed out that p, A 
and k of a real gas are vary with the temperature and density, e grows as # 1+7 with 7 ~ 0.5 and 
k increases like 9 q with q G [4.5, 5.5]. 

In the case that the data (po, 9q,uo, /, /i) are sufficiently regular and the initial density po is 
bounded away from zero, there exits a unique local strong solution to the problem (jl.ip - (jl.6p . 
and the solution exists globally in time provided that the initial data are small in some sense. 
For details, we refer the readers to papers O [T4"l \T7\ [20] and the references therein. When 
p, A, k are three constants, there have been some existence results on the strong solutions for the 
general case of nonnegative initial density. R. Salvi and I. Straskraba showed in [18] that if ft is 
a bounded domain, p{-) £ C 2 [0, 00), po S H 2 , uq S Hq n H 2 and the compatibility condition: 

— /idiv(Viio + Vud) — AVdivuo + Vp(po) = Po 9i f° r 9 £ L 2 , (1-13) 

is satisfied, then there exists a unique local strong solution (p, u) to the initial boundary value 
problem for the isentropic fluids. Independently of their work, H. J. Choe and H. Kim [7] proved 
a similar existence result when is either a bounded domain or the whole space, p = ap"'(a > 
0,7 > 1), p e L 1 n H 1 n W 1 ' 6 , u e Dl n D 2 and the condition (fLl~3|) is satisfied. Later, there 
are some results in [U [5] . 

Considering the physically important case that coefficients p, A and k are not constants, 
some results had obtained in one-dimension, see [HI EU [23j [25] etc.; Feireisl[l2] proved the 
existence of globally defined variational solutions to the compressible Navier-Stokes equations in 
W N with temperature-dependent viscosity; and some results had obtained for the incompressible 
Navier-Stokes equations, see (6[ El H5] etc.. In [6j, Cho-Kim study the case that (p,X,k) = 
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(p, A, K)(p, p9), and this method is fail to our case that (n, A,k) = (//, A, K)(p,9), since the 
density will be degenerate. In this paper, we shall combine the ideas in [5J |6l U [121 CES] to 
study strong solutions to the initial boundary value problem (jl.lj )- (jl.6p in a domain f2 C M 3 , 
with nonnegative initial density and (e,p, k, p, X) = (e,p, k, p, X)(p,9). The lack of simplicity 
of real flow (e,p, k, p, A) = (e,p, k, p, X)(p,9) makes the analysis significantly different from the 
ideal flow ( p = Rp9, e = C v 9, and C v ,k, p, X = constants). 

Throughout this paper, we use the following simplified notations for the standard homoge- 
neous and inhomogeneous Sobolev spaces: 

U = L r {Q), W k ' r = W k ' r {Q), H k = W k ' 2 , 

D k > r = {ve L} oc (n) : \\V k v\\ L r < oo}, D k = D k > 2 , 
Dl = {v G L 6 : \\Vv\\ L 2 < oo and v = on dQ}, H% = D\ n L 2 . 
It follows from the classical Sobolev embedding results that 

IHlis < CH^II^x, \\w\\ L oo < C\\w\\ w i, q , \\w\\ L oo < C\\w\\ D i nD i, q 

and 

\\iv\\l<=° < C\\w\\ D i nD 2, IMU°° < C\\w\\ L 6 nD i, q , 

provided that q > 3. Hereafter, we denote by C a generic positive constant depending only on 
the fixed constant q and the norms of (h, /). When Q is a bounded domain in M 3 with smooth 
boundary, we choose 9 C G C°°(R 3 ) satisfying 9 C = 9b near <9$7. When O is the half space or 
an exterior domain in M 3 with smooth boundary, we choose 9 C G C°°(1R 3 ) satisfying 9 C = 9b near 
dQ and 9 c \b c = 9°° where = {\x\ > R} and R is large enough. When O = M 3 , we choose 
9 C G C°°(R 3 ) satisfying 9 c \ B c r = 9°°. 

Our existence and uniqueness result can be stated as follows. 

Theorem 1.1. Let p°°,9b, 9°° > and 3 < q < 6 be fixed constants. Assume that (po, 9q, no, h, f) 
satisfy the regularity conditions 

Po > 0, po - p°° G L 6 n D 1 n D 1,q , (9 - 9 C , n ) G D\ n D 2 , (1.14) 

(h, f) G C([0, oo); L 2 ) n L 2 (0, oo; L q ) and (h t , f t ) G L 2 (0, oo; H' 1 ) (1.15) 
and the compatibility condition 

i 

-div(K V6» ) - Q(po,9 ,u ) = pfigi, - n Q ^ ^ 

L(p ,9 ,uo) +Vp{po,9 ) = p 2 g 2 , 

for some (51,52) G L 2 , where kq = k(pq,9o), Q(p,9,u) = ^^-\Vu + Vn T | 2 + X(p, 6>)(divn) 2 and 
L(p,9,u) = — div(u(p, #)(Vn + Vn T )) — V(A(p, #)divn). Then there exist a small time T* > 
and a unique strong solution (p,9,u) to the initial boundary value problem f 1. 1\) - (T1)\) such that 

p — p°° G C([0, T*]; L 6 H D 1 n D 1 ' 9 ), p t G C([0, T*]; L 2 Pi L 9 ), (1.17) 

(9-9 c ,u) G C([0,T*];Dj nD 2 )nL 2 (0,T*;D^), (1.18) 
(0 t ,«t) G L 2 (0,T,;D^) and (/M,p^n t ) G L°°([0, T*]; L 2 ), (1.19) 
l|p(*r)-po(OIU6 n0 ino^ + ll(^r)-^o(0,«(ir)-«o(0)llD n W^O, as i-0. (1.20) 
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Furthermore, we have the following blow-up criterion: IfT* is the maximal existence time of the 
strong solution (p,0,u) and T* < oo, then 

limsup(||p(t, •) - p°°\\ifinDinDi* + \\\fp~ e t{t, OIU 2 

t/T* 

+ j \\VO t {s, -)f L2 ds + \\9(t, •) - c |b oW ,« + IK*, -)bj) = oo- (1-21) 

Remark 1.1. The L 6 D -D 1 fl -D 1,9 -regularity of po — p°° seems inevitable to prove the local well- 
posedness in the framework of Sobolev spaces for any compressible fluid model in three dimen- 
sions, because the Sobolev embedding L 6 n D l ' q <^-> L°° holds only for q > 3 in M 3 . Moreover, we 
allow pq and p°° to vanish and so we may consider both interior vacuum and vacuum at infinity. 
Here, we only assume that po — p°° 6 L 6 fl D 1 n D ,q , which is weaker than the assumption in [6]. 

Remark 1.2. To simplify the presentation, we only give the proof of the case that Ob = 0°° = 
C = 0. 

Remark 1.3. As has been observed in [3j El [7J [18], the lack of a positive lower bound of 
po should be compensated with the compatibility condition (|1.16p . If (p,0,u) is a sufficiently 
smooth solution to (|l.l|) - (|1.3p . then letting t — ► in the equations (|1.2|) - (|1.3p . we easily derive a 
natural condition: there exists a pair ((71,52) of scalar and vector fields such that 

- div(K V<9 ) - Q(po, Go, uq) = Po9i and L(p , , u ) + Vp(p , ) = p g 2 in Q. (1.22) 

But it turns out that the weaker condition (11.160 is sufficient to prove the local existence 
and uniqueness of strong solution. Roughly speaking, the compatibility condition (|1.16p is 
equivalent to the L 2 -integrability of ^fpOt and ^fput at t = 0. Consider the elliptic system 
with variable coefficients, we obtain some regularity results in Section [5l Under these elliptic 
regularity results and the condition (11.160 . we could deduce that (Ot,ut) € L 2 (0,T; Dq) and 
(s/pOt, y/put) £ L°°(0,T; L 2 ) for a small time T > 0. Here, the compatibility condition (|1.16p 
is satisfied automatically for all initial data (po,@o, u o) with the regularity (11.140 whenever po is 
bounded away from zero. 

Remark 1.4. Consider the isentropic case and p, A being constants, B. Desjardins [10J prove the 
local existence of a weak solution (p, u) with a bounded nonnegative density to the periodic 
problem as long as 

sup (\\p(t)\\ Lao( j3) + ||Vu(i)|| L 2 (T 3)) < 00; 

0<t<T* 

Y. Cho, H. J. Choe and H. Kim [3] proved the local existence of a weak solution (p, u) as long as 

limswp(\\p(t)\\ H i nW i, q + \\u(t)\\ D i) < 00. 

t/T* u 

Using similar arguments in our proof, we can see that in the special case that the function 
ee = e# (p) is independent of the function 0, (jl.2ip may be replaced by 

limsup(||p(t, •) - p 00 \\ifinD 1 nD 1 « + \\( e ~ G c,u)(t,-)\\ D i nD i, 6 ) = 00. (1.23) 

t /T* u 

Furthermore, since the local existence time T* and all regularity estimates of the strong solu- 
tion in Theorem 11.11 depend only on || (51, <72)||z,2 and the norms of the initial data (po, Oo, uq, f, h). 
The proof of Theorem 1 1 . 1 1 also shows that the continuous dependence of the solution on the initial 
data holds at least for a small time interval. We may state the following result without a proof. 
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Theorem 1.2. For each i = 1,2, let (pi,9i,Ui) be the local strong solution to the initial boundary 
value problem hi. 1]) - [L6\) with the initial data (poi,0oi,uoi), which satisfy regularity conditions 
\1.14\ )- {1-15\) and the compatibility condition U.16\) with (51,52) = (gn,92i)- Assume also that 
for each i = 1,2, 

\\P0i ~ / 3 °°llL6nD 1 nD 1 .9 + \\(0Oi - 6c,U0i)\\Di(lD 2 + 1 1 (ffli) 92i) \\l 2 

+ W(h,f)\\L FL lnL*L% + W(htJt)\\ L 2 H -i < K 

andp = poi -P02 G X, (y/pbl{9 i - #02), y/m(uoi ~ U02)) G L 2 , where X = L 2 PiL%, if p°° = 0, 
X = L 2 , if p°° 7^ 0. Then there exist a small time T** and a positive constant Ck, such that 

sup (||p||i- + \\^/~0f L 2 + ||v^*ll|a) + / (WWdI + \Mdi) dt 
te[o,T.»] Jo v oy 

< c k (IIpoIIx + ll-^Poi^ollia + \\Vpo1u \\ 2 l2 ) , 
where ~p = pi — P2, = 6\ — 62 and u = u\ — Ui- 

The global result of the initial boundary value problem (ll.lj) - (11.6p with po ^ is very difficult 
to obtained. For a special case that the temperature vanish in the vacuum state, we can give a 
blow-up result in Section [6J 

In Section El we consider the following equations for a compressible fluid in W 1 x > 1): 

Pt + dw(pu) = 0, (1.24) 

(pu) t + div(pn (g> u) + Vp = div(p(Vu + Vn T )) + V(Adivn), (1.25) 

(pe) t + div(peu) + pdivu = ^| Vu + Vu T \ 2 + A(divu) 2 + div(«Ve), (1.26) 

where (k,p,X) = (k, p, A)(p, 9). In Section [6l we will consider only polytropic gas, so that the 
equations of the state of the gas is given by 

p = Rp8, e = c v 6, and p = Aex.p(S/c v )p' y , (1-27) 

where R > is the gas constant, A is a positive constant of absolute value, 7 > 1 is the ratio of 
specific heat, c v = is the specific heat at constant volume and S is the entropy. In Section 
[6l we only assume that 

/i(0,0)>0, nA(0,0) + 2^(0,0) = A > and k(s 1 ,s 2 ) > 0, (1.28) 

for all si, 82 > 0. The Navier-Stokes system is supplemented with the initial data 

{p,S,u)\ t=0 = (p ,S ,u )(x) G H k (R n ), k > [^] + 2. 

The blow-up of smooth solutions of compressible Euler equations (p = A = k = 0) has been 
studied by several mathematicians. T.C. Sideris [19] showed that the life span T of the C 1 
solution to the compressible Euler equations is finite when the initial data have compact support 
and the initial velocity is sufficiently large (super-sonic) in some region. When the initial density 
and velocity have compact supports, T. Makino, S. Ukai and S. Kawashima [16] shown the blow- 
up of smooth solutions to the compressible Euler equations in M 3 without external force and heat 
source. In a different way from [16] and [19], Z. Xin [22] showed the same blow-up result for the 
compressible Euler equations. Considered the case that p, A are constants and n = 0, Xin|22] 
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also showed similar results for the compressible Navier-Stokes equations for polytropic gas, when 
the initial density has compact support. In [3], authors extended Z. Xin's result [22] to the case 
of k > and fj,, A being constants. But their point of view in [U [22] cannot be applied for the 
problem in Section (6[ since the fact that 

/ [div(/x(Vu + Vu T )) + V(Adivu)] • xdx = 
it™ 

is strongly necessary in their argument, which seems hard to be obtained here. Here, we extend 
the Xin's blow-up result to the case that p,, A, k depend on the density and the temperature. 

Before stating our blow-up result, we introduce some notations. We denote by B r = B r (0) 
the ball in M. n of the radius r centered at the origin. We will use several physical quantities: 

m o = / podx (initial total mass), 

JR™ 



m \ = / po{x)xdx, 
^2 = / po{x)uo{x)dx (initial total momentum). 

JR™ 

We always assume that mo > 0. 

For the proof of the blow-up, we have only to prove the following theorem. 

Theorem 1.3. We assume that (p,\,n) satisfies \1.2S\) and T > 0. Suppose that (p,S,u) £ 
C 1 ([0,T];iI* : (R n )), k > [f]+2, is a solution to the Cauchy problem (L2^>-(T2^> with the initial 
data (po, Sq,uo), where the initial density po is compactly supported in a ball B ro . Then we have 

m r > mu + m^T > -m r , i = 1, . . . , n. 

Remark 1.5. If T* is the maximal existence time of the solution (p,S,u). Then since mo is 
strictly positive, this theorem implies that T* should be finite when m2 7^ 0. This theorem also 
shows the relationship between the size of support and the life span. For example, the range of 
life span can be extended as the support of the initial density become larger. Especially, we can 
expect the global existence of the smooth solution to the compressible Navier-Stokes equations 
in the case that the initial density is positive but has decay at infinity. 

Remark 1.6. As in [HE2], ° ur method need that the temperature vanish in the vacuum state, so 
we consider the entropy S and equation (|1.27j) which guarantee the temperature vanish in the 
vacuum state. 

The rest of this paper is organized as follows. First, based on some regularity results on 
elliptic system in Section [5j we obtain some a priori estimates for the linearized problem in 
Section [2j In Section [31 we will prove the existence results for the linearized problem. Then, 
using the usual iteration argument, we will prove Theorem 11.11 in Section [4j In Section [6l we will 
prove a blow-up result. 

2 A priori estimates for the linearized problem 

In this section, we consider the following linearized problem 

p t + div(H = 0, (2.1) 
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pu t + pv-Vu + L(a, 7], u) + Vp(p, 9) = pf, in (0, T) x O, 



(2.2) 



eg (a, rj) {pOt + pv ■ V0) — div(«(cr, rf)V9) + p(/?, 6*)divu 
= <5(cr, r/,u) + p/i + e p (a, rf)p 2 d\vv, (2.3) 

(p, 0,u)\ t =o = (po,0 o ,u ) in O, (2.4) 

(9, u) = (0, 0) on (0, T) x df2, (2.5) 

(p, e, u)(t, x) — » (p°°, 0, 0) as |x| — > oo, (2-6) 

where <t, 77, u are three known functions on (0, T) x Vt. Recall that Q(a, r/, v) = \ Vv+X7v T | 2 + 
A(a, 77)(divu) 2 and L(ct, 77, u) = — dw(p(a, r/)(Vu + Vm t )) — V(A(<7, rj)dwu). 

Throughout this section, we assume that the known data satisfy the following regularities 

Po > o, Po - P °° g L e n D 1 n D l,q , (e , u ) e d\ n D 2 , (2.7) 

(h, f) G C([0, oo); L 2 ) n L 2 (0, oo; L 9 ), (^, f t ) G L 2 (0, oo; H' 1 ), (2.8) 

(r, 1 ^c([o 1 r];Djnz) 2 )nL 2 (o 1 r;D 2 ''') ) fe.^e^to^;^), (2.9) 

eCdO.Tli^n^nD 1 ' 5 ), <r t G C([0,T];L 2 ni' ; ),ff% G L°°([0, T];L 2 ), (2.10) 

for three constants p 00 ,^ 00 > and q G (3,6). Then the existence of a unique strong solution 
(p,6,u) to the linearized problem (|2.ip - (|2.6p can be proved by standard methods at least for 
the case that po is bounded away from zero, i.e., po > 5. Using similar arguments in [U EJ [8], 
applying the classical method of characteristics, a semi-discrete Galerkin method and the usual 
domain expansion technique, we can obtain the following two lemmas. 

Lemma 2.1. Assume that po and vsatisfy the regularities 

Po>0, po-p°° GC andv G L°°([0,T}; n D 2 ) n L 2 (0,T; D 2 ' q ) 

for two constants p°° > and q G (3,6), where the space Co consists of continuous functions on 
f2 vanishing at infinity. Then there exists a unique weak solution p in p°° + C([0,T];Co) to the 
linear hyperbolic problem 112. 1\) . {2.J$ and 112. 6]) . Moreover the solution p can be represented by 



p(t,x) = p (U(0,t,x))exp I- J divv(s,U(s,t,x))ds^ , (2.11) 

where U G C([0,T] x [0, T] x £1) is the solution to the initial value problem 
f iU(t,s,x) = v(t,U(t,s,x)), te[0,T] _ 

\ U(s,s,x) = x, se[o,T],xen. K ' ' 

Assume in addition that po — p°° G L e D D 1 n D ,q . Then, we also have 

p — p°° G C([0, T]; L 6 n D 1 D D 1,q ) and p t G L°°([0, T];L 2 D L q ), (2.13) 

and 



/ ?0 °llL 6 nD 1 nZ) 1 ><? < (^\\Po ~ P°°\\L f 'nD 1 nD 1 ''i + Cp 00 \\v(s,-)\\ D 2 nD 2, q ds 

xexp^C J \\Vv(s, ^Ww^nD^ds^j (2.14) 
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Lemma 2.2. Assume in addition to \2. 7\ )- [2~Id]) that po > 5 in £1 for a constant 5 > 0. Then 
there exists a unique strong solution (p,0,u) to the initial boundary value problem \2. l\) - (2Zh)) in 
[0, T] such that 

p — p°° £ C([0, T]; L 6 n D 1 f] D 1,q ), p t £ C([0, T]; L 2 n L q ), (2.15) 

(0, u) G C([0, T];Dq fl D 2 ) n L 2 (0, T; D 2 ' 9 ), (2.16) 

(0 t ,«t) GL 2 (0,T;^) nC([0, T];L 2 ), (2.17) 

p > 5 on [0, T] x H /or a constnat 5>0, ( 2 -18) 

||p(i, •) - poQIbaDW^ + •) " W,u(t, •) - uo(-))||ujn£»!> °> as t 0. (2.19) 

Assume that the initial data (po,9o,uo) satisfy the hypotheses of Lemma 12.21 and the com- 
patibility condition 



(2.20) 



[ -div( K (o-(O),r ? (O))V0 o )-Q(o-(O), i ? ? (O), V (O)) = p> 9l , 
\ L(a{0),7 ] (0),u )+Vp(p ,e ) =p\g 2 in Q. 

Let us choose a fixed cq so that 

co > 1 + P°° + \\po ~ P°°\\l^d^d^ + ll(#o, u ) || D in£»2 + ll(0l>0a)ll£a, 

and 

c > + ||cr(0, -) -(T^IUenBinDi,, + ||(r?(0, •), w(0, -))|| D i nD2 , 

where 3 < g < 6 and o -00 ,/? 00 > are three constants. Moreover, assume that o-,rj,v satisfy the 
regularity stated in Lemma 12.21 and 

sup \\v(t, 0|U + / (\\v t (t, -)\\ 2 + \\v(t, .)|H, a ,,) dt < ci, (2.21) 
te[o,T] Jo v y 

sup ||«(t,0||j^ < c 2 , (2.22) 
te[o,T] 

sup (|| (V^Vt)(t,-)\\ L 2 + \\v(t,')\\m)+ I \\Vt(t,-)f D idt<c 3 , (2.23) 
te[o,T] jo 

SUp (||f7 t (t, OlbnW + ||<7(t, " <7°°|li«niJW.«) ^ C 4> ( 2 - 24 ) 

te[o,T] 

for some fixed constants c\, C2, C3 and C4 such that 

1 < c < c\ < c 2 , 1 < c < c 3 ,c 4 . 

Then we could derive some a priori estimates for the solution (p, 9, u) which are independent of 
5. 

I. Estimates for the density and the coefficients (//,, A, K, eg, e p ). 
To estimate the density p, we first recall from (|2.14|) that 

IIp(V) - P^Wl^d^d^i <Cc exp^C^ HV^OIIifinw^sJ 
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for t £ [0, T]. Then, using Holder's inequality, we obtain 



\\^ v \\H 1 nw 1 'ids < t2 
and conclude that 



\\Vv\\ 2 HlnWl , q ds 



< C(c 2 t+(c 2 t)2), 



l|p(*,-)-P oo ||r6nDW.«<Cc , (2-25) 
for < t < min(r, To), where To = c^ 1 < 1. Moreover, it follows from (|2.1ip and (|2.25j) that 

C~ x 8 < p(t,x) < Cc for < t < min(T,To),x G O. (2.26) 

From (|2.2ip - (|2,23p . using the Gagliardo-Nirenberg inequality, we have 

\\v(t) ~ v(0)\\dI < C f htWulds < C (J* IMI^s) 2 tk < C&K (2.27) 



6 — q 3q — 6 Sq — 6 
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\m-vm\D^<CHt)- V (0)\\^ o \W)-v(0)\\ D l q <Cc^ t—, (2.28) 



\\v(t) - v(0)\\ D i <CJ^ \\v t \\ D ids < Celts, (2.29) 
||t7(t) - vCO)!!^ < C]|«(t) - «(0)l|^"||t;(t) - vCO)!!^ 2 < Cc^c^t^, (2.30) 

Aq 5q-6 4q' gg^lg 

for (/ = max(<7,4). Let Ti = min(To, CgCg , Cq~ 9 c 3 6-9 , CqC^ 1 , Cq~ 9 c^ 1 c 2 6-9 ), we have 

[I^OIIuJniJM + II^OIIjjinijW <Coo, V0<t<min(T,Ti). (2.31) 
Using Sobolev's embedding theorem, we have 

|| (77, «)(*,-) || L oc <Cc , V0<i<mm(T,T!). 
From ([S3|> . (I2T25]) and (|2~3T|) . we obtain 

+ WptWvm* < Aocl (2.32) 

where < t < min(T, Ti). Here, the constant Aq > 1 is independent of any one of Cj, i = 0, . . . , 4. 
Using the interpolation inequality, we have 

\\p(t, •) - po(-)IU« < C||p(t, •) - PoOlli IM*, •) " Po(0ll!~ < 

where < t < min(T, Ti). 

Choosing c 4 = Aq(%, from (fTTOj) - (fl~TT|) . (l2T23ll - (l2T24l and (I23TD . we have, for any < t < 
min(T,Ti), 

\\VfM(a(t,-),ri(t,-))\\ L2nLq 

< \\n P (a(t, -), V (t, -))Va(t, -Jll^nw + IM*(*> O^fo 0)V»#, Oll^nM 

< C(co) (||V<r(t, -)\\ L 2 nLq + ||Vr?(i, OlbnL<0 < C(c ), 

and 

||/^(a(t,.),»7(t, 0)11x3 
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< \\fji p (a(t, -),ri(t, -))a t (t, -)\\ L 3 + \\p e (a(t, •),#, ■))%(*, 

< C(cq) -)||z,3 + lk^77 t (t, OlUa) 

< C(c )\\a t (t, -)\\ L3 + C(c )\\a^ Vt (t, -)\\l 2 \Mt, -)|||a 

< C(c )(c 3 + c|||r7i(t,-)ll|i)- 
Similarly, we obtain 

sup (||(//,A,K,e e ,e p )(a(V),^(V))llcr/3 + II (Po, Pi, P2) -))IU°°) ^ c 5> ( 2 -33) 

0<t<min(T,Ti) 

sup ||V(/i, A, k, Bo, ep)((j(t, •), r](t, •))HL 2 nL<? — c 5) (2.34) 

0<t<min(T,Ti) 

sup \\dt(n,X,K,ee,ep)(a(t,-),r](t,-))\\ L 3 < c 5 ( c 3 + cf ||»ft(t, J . (2.35) 

0<i<min(T,Ti) V °/ 

for a constant C5 > 1 which depend only on Co, and /?= 1 - J E (0,1). 
II. Estimates for the internal energy and the pressure. 

Taking the operator dt to the equation (12, 3j) . we obtain 



eepfltt - div(«;V6' i ) - div(«iV0) + (pdivu) t 
= Q(a, r], v) t + {ph) t - (egpv ■ V9) t - {e e p) t 9 t + (e p p 2 divv) t . 

Then multiplying this by 9t, integrating over f2 and using (|2.ip . we have 

\~[ t \ e eP^dx + j n\V9 t \ 2 dx + j (pdivv) t 9 t dx - J 9 t di\(n t V9)dx 
= J (~(e p) t Ot +0t{-(e e pv ■V9) t + Q(a,ri,v) t + (ph) t + (e p p 2 divv) t )) dx 



and 



\jt j e eP 9 t dx + K ° / Mt\ 2 dx 



< c /(i^iiv^ii^i+pIv^ii^i + i^iiv^iiv^i + h^ii^i 



2 

+\pv ■ V(e e 9 2 )\ + \(e eP v ■ V9) t \\9 t \ + {p+ |A|)|Vu||Vi*||0 t | + (M + |At|)|V«| 2 |0*| 

11 

+\pt\\h\\9 t \ + \h t \p\9 t \ + \9 t \\(e pP 2 divv) t \)dx :=J2 T r ( 2 - 36 ) 

i=i 

Using (|2~2T|) . (l2T25]) - (p^6l) and (p32l) - (l235| ). applying the Cauchy- Schwartz, Gagliardo-Nirenberg, 
Holder and Sobolev's inequalities, we can estimate each term Ij(t) for < t < min(T, Ti) as fol- 
lows: 



h < C (\p e 9t\ + \p P pt\)\Vv\\9 t \dx 

< C J Pl (l + 9)p^\9 t \ 2 \Vv\dx + C J p 2 (l + 2 )|p t ||Vt;||0 t |dc 

< Cc 5 (||fl|U«||V«|| L 6 + \\Vv\\ LS ) \\Vp9 t \\ L 49 t \\ Dh 
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+Cc 5 (||(9||i6||Vt;|| L 6 + ||V V || L2 )lhl|L3||e t bi 

< c44 (i + ||^|| 4 D(1) ) (i + ii^lli.) + ^liv^||| 2 , 

h < c J po(i + e 2 )p^\w t \\e t \dx 

< C C 5||^IIl 2 ||V^|| l2 (i + ||0||1oc) 

< iiv^ni, + c4\\^pe t f L2 (i + \\e\\* Dh + \\9\\% 2 ), 

h < c\\ Kt \\ L s\\ve\\ L 4v9 t \\ L , < c4(4 + c 3 \\ Vt \\ Dh )\\9\\l 2 + ^l|w t ||| 2 
< c44\\ef D2 + ^11^11^ +c c fc|||^|| 4 D2 + ^||v^||| 2 , 



< c44\\^pe t f L2 + ^11^1, 

h < C {\\d t ee\\ L 4p\\ L ~\\v\\L~ + \\eo\\L«>\\pt\\ifi\\v\\L<*>) \WO\\lAMdI 

H-C7||e < »||ioo|| >£ ,||2oo |IV^»*IU=*||vt|| J D 5 l|V0||i3 H- ||e tf |U~||^||ioo||v||x.oo|| > /^*IU=»ll^*|U 2 

< c4(4 + C3i|r ? ,||^) c 4 ||v^|ii 2 + c44\\^- P 9 t \\l 2 + ^||V^||| 2 
^v t \\l h +Cco4\\ve\\ L 4e\\ D 4^pe t \\l 2 , 

h < Cc 5 ||Vt;|| L 3||V^|| L2 ||^|| L 6 < C44\\Vv t \\ 2 L2 + |j||V0 t ||| 2 , 

/ 8 < C(ll^llL3 + l|A t ||L3)||(V^) 2 || L2 ||^|| L6 

< Cc 5 (c 3 + c\\\r ]t \\l h )4\\V9 t \\ L , 

< C44(4 + c 3 \\r, t \\ Dh ) + |j||W t ||| 2 , 

h < c\\p t \\ L -4h\\ L2 \\9 t \\ L , < c4\\h\\ L 4V9 t \\ L , < c4\\h\\ 2 L2 + ^\m\\l 2 , 
h < c\\h4 H -4 P 94 H i 



h < cii^iUsii^ll^lblllocll^bi 

< cc 4(4 + c 3 Ki| D i)||^ t |li 2 + ^liv^||| 2 , 




< 



C (\\p\\U\v\\ L ~\\Ve e \\ L 3 + ||p||ioc||«||L<»||e 9 || L =o) \\y/p6t\\v\\et\\ Dh 
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1 



< c||/Hllif-i(||p|lloo||^0t|| L2 + \\ P \\L^e t \\ L , + ||Vp|| L3 ||v^|| L 2) 

< c4\\h t \\ 2 H ^ + \\ p h t \\l 2 + ^\\ve t \\l 2 , 



and 



111 < ll^lUallpll^HV^II^II^II^ + ||eJ Loo ||p|MH| L 3||V<;|| L2 |^|| D i 
+l|eplk-||p||I oc ||V^|U 2 ||Vp^ll^ 
< Cc|(cI + C3|| % || B pc« + ^||^||^ + ||V Ut ||| 2 + ^ci||^dl!- 
Substituting these estimates into (|2.36|) . we have 

j t \\VwVt\\h + ^e t \\l. 



< c44(i + ||w||£ 2 + ||^|| 4 D2 )(i + \\ P ^ t f L2 ) + c44\\e\\l 2 + 1 



+C44\\9\\% 2 +Ccgci||v t ||J, 5 +c44(4 + c 3 \\ Vt \\ Dl )(\\^- P e t \\l 2 + \\V9\\ 2 L2 ) 



+c c y 5 (4 + c 3 \\rj t \\ Dh ) + c4\\h\\i 2 + c4\\h t \\^ 

< Cc 8 2 44(l + ||W|| 4 i2 + \\9\\% 2 )(1 + ||/M|| 4 L2 ) + CC3 + CcgcjjHIJg 

+c4 2 44 + c4\\h\\ 2 L2 + c4\\h t \\ 2 H ^, 

and then integrating this over (r, i), using the estimates (|2.8p . (|2.2ip and (|2.23p . we have 

\\V^pO t (t)\\ 2 L2 + J' \\V9 t \\ 2 L2 ds 
< C\\(^9 t )(r)\\ 2 L2 + C4c lC 2 + Cc l2 44t 

+C444 f\l + ||V0||i 2 + ||0|| 4 D2 )(1 + \\ P l i9 t \\i 2 )ds, (2.37) 

for < r < i < min(T, Ti). To estimate limsup r ^ \\y/^ep9t{T) <[ \ 2 L2 ^ we observe from the equation 
(El that 



? 2 dx < C j (p\h\ 2 + e e p\v\ 2 \V9\ 2 + p^p^Vvl 2 

+p- 1 (div(KV9) + Q(<t, r?, v)) 2 + e 2 p 3 \ Vv\ 2 ) dx 
and using the estimates (12.20P and (I2.33h . we obtain 
limsup \\^/eep9t(T)\\ 2 L2 



Hence, letting r — » in (|2.37p . we deduce that 

\\pt0t{t)\\h+ fwwtWhd* 

Jo 
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< Cc 8 2 44 [ (l + \\V0\\l2 + P\\%z)(l + \\P^t\\h)ds + C4cx4, (2.38) 
Jo 

for < t < min(T, T2), where T2 = min(Ti, Cq 8 Cg 3 c 5 -2 ). On the other hand, since 

|V0||| 2 =2 J ve ■ ve t dx < 2\\ve\\ L 2\\v9 t \\ L 2, 



— \ 
(ft 



it follows that 



||V0(i)||£ 2 < c\\ve f L2 +c ! \\ve t f L2 ds, 

J 

for < t < min(T, T2). Combining this and (|2,38p . we obtain that 



\\V0(t)\\h + \\P*0t(t)\\b+ / IWtWfrds 

Jo 

< Cc s 2 44 [ (l + ||V0||t 2 + ||0|^ 2 )(l + ||p^t||| 2 )d S + C^ci^, (2.39) 
J 

for < t < min(T,T 2 ). 

To estimate ||0||£>2, we use the following elliptic regularity results: if (/, g) G Dq n JD ' 3 is a 
solution to the elliptic system 

-div(«V/) = F and - div(/i(V# + V# T )) - V(Adiv 5 ) = G, mQ 

with (F, G) e L 2 HL q . From (|2"33>(j2TM]) . we have 

3 

sup \\(K,fi,X)(a(t,-),r](t,-))\\ C (3 nD i, q < c 5 , j3 = l--. 

0<t<min(T,T 2 ) Q 

From Lemma 15.51 and 15.71 we have (/, g) G D 2 satisfying 



\o*(n) < c 6 [\\F\\ L 2 (n) + ||V/|| i2(c) + ||VK||2- ( 3 n) ||/|| D i (Q) J , (2.40) 

llsllz^o) < c 6 (l|G|| L2(n) + ||V 5 || L2(n) + |||Vm| + iVAHI^II^IIoi^^ , (2.41) 
and from Lemma 15.61 and 15.81 we have (f,g) G D 2 ' q satisfying 



2q-6 



D 2, Hn) < c 6 [\\F\\ Lm + ||V/|U 9( n) + \\^4lT { n) ll/ll^(n) ) , (2-42) 

\\9\\i»«{n) < c 6 (||G|| L9(n) + ||V<7|| i9(n) + |||V/x| + IVAIH^ \\g\\ D i m ) , (2-43) 

where the constant cq = cq(co) > 1 depends only on cq. It should be noted that the estimates 
(|2,40p - (|2,43p hold for both bounded and unbounded domains. Applying the elliptic regularity 
result (12.4Q[) to the equation — div(«;Ve) = F in Q, where F = p(h — egOt — egv ■ V9) — pdivv + 
Q(a,r),v) + e p p 2 divv, using Holder and Sobolev's inequalities, we have 

||0|| D 2 < Cc e (\\ph\\ L 2 + IMIlHIp^IIl 2 + ||e e ||i,oo||pi; • V0|| L 2 + ||pdivu|| L2 

+ \\Q{a,ri,v)\\ L 2 + ||V0|| L 2 + || Vk||£~ 3 q) ||0|| D i(n) + ||ejz,°° \\p\\ 2 La o || Vu|| L 2) 
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< Cc 3 c 5 c 6 (l + \\p*6 t \\v + \\V9\\ L 2) + Cc 6 cf 3 \\V6\\ L 2 
+Cc 6 \\po^p\\l~\\(1 + 2 )Vv\\ L 2 

< C4cf~ 3 c 6 (l + \\ph t \\ L 2 + \\V9\\ 2 L2 ). (2.44) 
Combining this and (12.39p . we deduce that 

||v^(t)||| 2 + || P ^ t (t)||| 2 + f \\ve t \\l2ds 

Jo 

< Ccf4cf~ 3 4 A(l + ||V0||| 2 + \\ph t \\^)(l + \\ph t \\t 2 )ds + Cctc.cl (2.45) 
J o 

for < t < min(T, T2). Now we define a function H{t) by 

ff(i) = i + ll P ^ll! 2 + llv0|i! 2 . 

Then it follows from (pT35]) that 

H{t) < C4a4 + Ccf4cf 3 4 / H 6 (s)ds for < t < min(T,T 2 ). 

Jo 



Solving this integral inequality, we easily derive 



21a. j 

H(t) < c4ci4i l - c4 5 4 c r' i 4 t y* for a11 sma11 * > 0- 

llq_ 

Therefore, taking T3 = min(T2, (2C4 5 4 C 5~ 3 c i) _1 )> we conclude that 

\\o{t)f D} + wMmh + T iietWiii,!* < a.4^1 (2.46) 

and 

\\e(t)\\ D 2<A l 4cf~ 3 c 6 , (2.47) 

for < t < min(T, T3). Here, the constant A± > 1 is independent of any one of Cj, i = 0, . . . , 6. 
From (|2.46p - (|2.47p . using the Gagliardo-Nirenberg inequality, we have 



t / rt 



\\0(t) -0 o \\ D i <C J \\e t \\ D ids < C U \\9 t \\ 2 DlQ d S j t~2 < C4cfc 5 t2 < Ceo, (2.48) 



1 _ n i 1 



\\e(t)-e \\ D i, q < c\\e(t) - e \\ 6 *! \\e(t) - e \\}/ 



2q 



6-g 



< C4cf 3 c 6 t— < Cc , (2.49) 

_ 32q 8q 2 4q 

for < t < min(r,T 4 ), where T 4 = min(T 3 ,c 2 6-9 c 5 (^ 3 >( 6 -^c 6 6-9 ). Therefor, from (fl~71) -(T01). 
(|2.46|) and (|2.48p - (|2.49p . using Holder and Sobolev's inequalities, we obtain 



sup [\m,-)\\ Dhnm , g + \\e(t,-)\\ L ~ )<Cco, (2.50) 

0<t<min(T,T 4 ) V 7 

sup \\p(p{t, •), 0(t, <C4c 5 , (2.51) 

0<t<min(T,T 4 ) 
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sup \\v P (p(t r ),e(t,-))\\ L 2 nLq 

0<t<min(T,T 4 ) 

< sup (||PpV^|| L 2 nL9 + \\peV0\\ L 2 nLg ) 
0<i<min(T,T 4 ) 

< sup {\\P2(l + 2 )Vp\\ L 2 nLq + \\ Pl ^p(l + 9)Ve\\ L 2 nLq ) 

0<i<min(T,T 4 ) 

< Cc 3 c 5 , (2.52) 



and 



sup \\d t p{p{t, -),9{t, -))\\ L 2 < sup (\\PpPth 2 + WPoOth 2 ) 

0<i<min(T,T 4 ) 0<i<min(T,T 4 ) 

< sup (\\p 2 {l + 8 2 )p t \\ L 2 + \\ Pl J-p{l + e)e t \\ L 2) <Cc\cl (2.53) 

0<t<min(T,T 4 ) 

Hence by virtue of (12^21) . from (USD, (I235D . (12311) and (l2T33|) - (p3iD . using Holder's inequal- 
ity, we deduce that 

f* \\V0(8)f D1 , q d8 

Jo 

< Cc\ [ (\\ P h\\l q + lleelli. (\\pO t HI, + • V0||| s ) + ||pdm;|li« 

•/ o 

5g-6 

+\\Q(<r,V,v)\\h + l|e P ||i=o||p||ioo||V«||i, + ||V0||£, + || Vk|| l 7 ( 3 c) ||0|&i (n) )ds 

< AcolWll* + ^([Ip^tH^ + \\d t f D r) + cgc§ + ciUV^II^M + Clcf^]ds 
JO 

< Ccfacjcl (2.54) 

for < t < min(T,T 4 ). 

III. Estimates for the velocity. 

To derive estimates for the velocity u, we differentiate the equation (|2.2p with respect to t, 
multiply it by ut and integrate it over f2, using the equation (|2.ip . then derive 

™ / P\ Ut \ 2dx + y (/"[ Vn * + Vu t\ Vu t + X(divu t ) 2 )dx 

= j u t - [-Vpt + (pf)t ~ Ptv ■ Vu - pv t ■ Vu + div(pv)u t 

+div(p t (Vu + Vu T ) + V(X t divu))]dx. (2.55) 

Using the technique in [12], we can estimate the second term in the left hand side of (|2.55p as 
follows: 

Tl . T7„,. i \{A\w,.\2\ 



J (p[Vu t + VuJ } : Vu t + \(d\wu t Y)dx 

> J (p[Vu t + VuJ] : Vut - ^p(divu t f)dx 

^|Vn 4 + VuJ — -d\vutl\ 2 dx 

U° f 2 

> — \Vu t + VuJ d\vu t I\ 2 dx 

2 J 3 



15 



/ 

J (p,°\Vu t \ 2 + ^°(dwu t ) 2 )dx > J fi°\Vu t \ 2 dx. (2.56) 



(fi°[Vu t + Vu t T ] : Vut - \n Q (d\vu t ) 2 )dx 



Using Holder, Young and Sobolev's inequalities, we easily deduce that 
|||p^|l! 2 + M°l|V^||l 2 

< c(\\ Pt \\ 2 L2 + IMMIp^HI* + (i + llPlli- + \\Vp\\b)\\ft\\ 2 H-i + llftlli.il/lli, 

+\\Pt\\h \Hlao || Vuf L2 + \\p\\ LaB \\v\Hco \\p$ut Ilia) + Cci\\p\\ LO o \\Vu\\ l2 \\Vu\\ d i 
+cr 1 ||V^|li 2 ||^u i ||i 2 +C(||^|||3 + ||A t ||i 3 )||n|| 2 D2 . 

Then it follows from the estimates (|2jgj> . (I235D . (l2T3TD - (l2T32ll . (12351) and (l233|) that 

^llp'Wtllia + IIVutHi, 

< cc\4 + ccgn/tii^: + (c4 + c^iiv^ni^iip^tiil. 

+C(cg + 44)Mh + Cc 2 (c 2 + c 3 ||^(t, OIlDplkllD- 
for < i < min(T, T4). Hence, using the facts that 

/9(r)-5 [L(a(r), t/(t), «(r)) + Vp(p(r), 0(r))] -» 52 in L 2 as r -> 

and 

I2,i < C||«o|loi + C / l|Vn t || 2 2 ds for < t < min(T,T 4 ), 
J 



we easily obtain 



1 /"* 

m*)llm + IIP^II^+ / HVutlliarfs 
Jo 



< Ccl + C [ (4c 2 + c r 1 ||V Wt || 2 , 2 )(||n|| 2 1 + ||p^ t ||| 2 )d S 
Jo 



+cc 2 j (4 + c 3 \\r, t (t, OH^ONI;^, (2.57) 

for < t < min(r,T 4 ). 

To estimate ||u||o2, we apply the elliptic regularity result (|2.4ip to the equation (|2.2p . Then, 
it follows from IfEBjl . (12351) . (l23Tl) - ([232j) and (123511 that 

Nil* < Gc 6 (||p/|| L 2 + ||pM 4 || L2 + ||pt;-Vu|| L2 + ||Vp|| L 2 
+ ||Vn|| L2 + |||V^| + |VA||||7l|V U || L2 ) 

< CcbOIpHlcc ||/||^ + llplll^Hpl^lU. + \\p\\ L ^\\v\\ L ^\\Vu\\ L2 
+^c 5 + ||u|| D i+c|^||Vu|| i2 ] 

1 ! 1 

< Cc 6 [c + 4\\p2u t \\ L 2 + 4c 5 + (4 + cf 3 )\\u\\ D i] 

< C4c 5 c 6 + C4c§ = *c e (\\u\\ D i + H^utllia), (2.58) 
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for < t < min(T, T 4 ). Therefore, substituting this into (|2.57j) . we obtain 



u(t)\\l h + \\p*Mb+ L \\Vut\\hds 







Jo 
and using Gronwall's inequality, we obtain that 

Ht)\\h n + \\p*Mt)\\h + f \\n{s)f Dl ds < c4eMC + Cctclcf~ 3 4r*} < c4, (2.59) 

Jo 

and 

\\u(t)\\ D 2 < Cc 5 cf~ 3 cs, (2.60) 

for < t < min(T,T 4 ). 

Moreover, it follows from (p03|) that if < t < min(T,T 4 ), then 

q 5q-(j 

\\u\\ D 2, q < Cc 6 (c \\f\\ Lq + co(\\p*ut\\ L 2 + \\u t \\ D i) + c 7 Q cf 3 c 6 + 4 q - 6 \\Vu\\ L 2) 
< Cce(co\\f\\ L , + 4 + c \\u t \\ Dh + c 7 cf~ 3 c 6 + cf^cf ), 



and thus 



'o 

for < t < min(r,T 4 ). 

Combining (|2.59p - (|2.6ip . we finally conclude that 



/ \\u{s)\\ 2 D2 , q ds<Cclcl (2.61) 
Jo 



Dl + \\P^Mt)\\^ + I (||«t(«)||?,i + \Hs)f D2tq ) ds < A 2 c 7 cl (2.62) 

and 



o 



\\u{t,-)\\ D 2 < A 2 c 5 cf 3 c 6 , (2.63) 

for < t < min(T, T 4 ). Here, the constant A 2 > 1 is independent of any one of Cj, i = 0, . . . , 6. 
IV. Conclusion. 

Let us define c±, c 2 and c 3 by 

q 2q 

Cl =A 2 c 7 cl, c 2 = A 2 c 5 cf 3 c 6 and c 3 = A lC 8 2 cf 3 c 6 . (2.64) 

Then choosing any such that < < min(T,T 4 ), we conclude from (12321) . fliOBI) - (12^47]) . 
(12341) and (J2j62])-([2j)3]) that 

sup (\\p(t, •) - p°°\\ 2 L e nD x nD i, q + \\pt(t, Oll^nw) < c 4 , (2.65) 
0<t<T„„ 

sup (\\(^p6 t )(t,-)\\ L2 + \\e(t,-)\\ D 2)+ [ T *** \\6 t (t,-)\\ 2 D} dt<c 3 , (2.66) 



te[o,T,„] jo 
sup (||0(t,-)|U + \\Vpu t (t,-)\\ L2 ) + f T *" \\e(t,-)f D2 , q dt < Cc 7 , (2.67) 

0<t<T,„ V u / Jo 
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sup \\u{t, OH ni + / (\\u t (s, Oll^i + ||«(a, OH!**) < ci, (2.68) 



0<t<T«** 

sup ||u(t, OIId 2 ^ 5 ^ c 2 5 (2.69) 

0<t<T»»» 

where C7 = CqCic^Cq. Using similar arguments in (|2.27p - (|2.29[) . we have 

\\(0(t, - e (-),u(t, O - n o (0)|| D i < CcfP, (2.70) 

and 

l|p(V)-po(-)llr« <C<$tf, (2.71) 

for < t < T«*. Hence, it deserves to emphasize that c±, . . . ,cy and T*** depend only on cq, but 
not on the lower bound 5 of the initial density pq. 

3 Existence result for the linearized problem 

Now we can prove the key lemma to prove Theorem II. H 

Lemma 3.1. Assume in addition to (2. r /)- (2~Tu^) that the initial data (po,9o,uo) satisfy the 
compatibility condition 



1 1 
- div(K V6» ) - Q(po,6q,u ) = Pq9i and L(p ,6 ,u ) + Vp(p ,6 ) = p^g 2 , (3.1) 

for some (51,(72) G L 2 . Assume further that 

M,v)(0) = ( Po , 0o, u ), a 00 =p°°, 



sup ||«(* s -)]lx»5 + / (\Mt r )\\ 2 D} + \\v(t,-)\\ 2 D2 , q )dt< Cl 

t£[0,T t „] ° JO V ' 

sup \\v(t, -)\ B 1 < C 2 , 

te[o,T***] 

sup (\\(y/ffT k )(t,-)\\v + \\ri(t,-)\\ D 2)+ f \\vt{t,-)\\ 2 D idt<c 3 , 

te[0,T»„] JO 

sup (||cr(i, - P° ||LBnD 1 nD 1 .9 + \Wt(t, Oll^niO < c 4 , 
te[o,T»»»] 

/or i/ie positive constants c\, 02, c%, C4 and T^, depend only on cq, chosen as before, where 

c = 2 + p°° + ||/3 - P^IUonDino*.* + IK^o, w )|| D i nD 2 + \\(gi, 92)\\ 2 L 2. 

Then there exists a unique strong solution (p, 8, u) to the linearized problem 12.1]) - [276]) in [0, T* 
satisfying the estimates $2.65]) - [2.71 ) as well as the regularity 

p-p°° G C([0, T„»]; L 6 HD'n D 1 ' 9 ), p t £ C([0, T„*]; L 2 n L 9 ), 

(0, u) g c([o, r***]; Z)q n d 2 ) n l 2 (o, t«* ; d 2 - 9 ), 

(0 t ,«t) G L 2 (0,r^ ;J D^) and (p^p^t) G L°°([0, T***]; L 2 ), 
||p(V) -po(OIUe nJ DinOi.9 + II (#(*>•) - ^o(-)> «(*>•) - n o(0)llDinD2 0, as t -> 0. 
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Proof. We define p$ = po + 5 for each 6 G (0, 1). Then from the compatibility condition (|3.ip . 
we derive 

f -div(K( Po ,9 )V6 ) - Q(p o ,0 Q ,u Q ) = (p s )hi, 

\ L(po,e ,u ) + v P (pie ) = 

where 

8 fpoV* , 5 fpcA 1 , V(p(p S o,Oq)-p(po,0o)) 

9i = [ — ) 9i and ff2 = — g 2 + — ! . 

\PoJ \PoJ (pg)5 

Since p p ,pe £ C a ([0, oo), L^ c ([0, oo))), a G (^il)> w e obtain that for all small 5 > 0, 

co > 1 + {p°° + 5) + \\p s - (p™ + 6)\\ L e nD i nD i, q + ||(0o, n )|| D i nD2 + \\(gl,d)\\h- 

Hence, from previous results for positive initial densities, we deduce that corresponding to the 
initial data (pq,0q,uq) with small 5 > 0, there exists a unique strong solution (p s ,0 s ,u s ) of the 
linearized equations (|2.1j) - (|2.6p satisfying the local estimates (12.65|) -( |2.69p . From this uniform 
estimates on 5, we conclude that there is a subsequence of the solutions (p s ,0 s ,u s ) which con- 
verges to a limit (p, 9, u) in an obvious weak or weak—* sense. It is easy to show that (/?, 9, u) 
is a weak solution to the linearized equations (|2.ip - (|2.6p in [0,T***]. Finally, thanks to the lower 
semi-continuity of various norms, we find that (p,9,u) also satisfies the estimates (I2.65p - (l2.7ip . 
This proves the existence of a strong solution (p, 9, u) with the regularity 

p — p°° G L°°([0, T***]; L 6 H -D 1 H D 1,q ), p t G L°°([0, T„*]; L 2 n L q ), (3.2) 

(0,u) e L°°{[0,T^}; D x D D 2 ) D L 2 (0,T^; D 2 ' q ), (3.3) 

(p t ,ut) G L 2 {0,T^;D*) and (ph t ,p^u t ) G L°°(0, T***; L 2 ). (3.4) 

Now we prove the uniqueness of the solution to this regularity class. Let (pi,9i,ui) and 
(P2> 62,1*2) be two solutions to the problem f)2. lj) - f)2.6j) satisfying the regularity (|3.2p - (|3.4p . and 
we denote 

P = pi - p 2 , 9 = 01-02, U = Ui- U 2 . 

First, since ~p G L°°(0, T***; L°°) is a solution to the linear transport equation p t + divQcw) = 0, 
it follows from the uniqueness result by DiPerna-Lions in [TT] that p = 0, that is, pi = p% := p 
in [0, T***] x fi. Next, to show that = in [0, T***] x f2, we multiply the both sides of 

e e {p0 t + pv • V0) - div(«V0) + dwv{p(p, 0i) - p{p, 2 )) = (3.5) 

by and integrate over (0, t) x f2. Then, since pt + div(pv) = and 0(0) = in S7, we deduce at 
least formally that 



e e p\9\ 2 dx + J J n\V0\ 2 dxds 



2 

1 

2 
.0 







—2 —2 — 

d t eep0 +Veepv0 — 0dvvv(p{p, 0i) — p(p, 2 )) 



dxds 



U ft f ft 

- Tj J\VWdxds + Cj [\\d t ee\\ 2 L3 + \\Ve e \\ 2 L3 \\v\\ 2 Laa 

+[|V«[|i3[|pi||i<»(l + ||(ei,e2)||i»)] \\Vp0\\hds (3.6) 

and applying Gronwall's inequality, we also conclude that = in (0, t) x VL. But this argument 
is somewhat formal since it is not obvious that p^0 G L°°(0, T***; L 2 ) for the case of unbounded 
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domains. Hence we have to justify this formal argument by deriving f|3.6|) rigorously. For this 
purpose, we assume that Q is an unbounded domain, let U r = Pi B r , r > 1, and define 
9 r G L°°(0,T^;H^{n r )) by 

9 r (t,x) = d(t,x)4> r (x) for (i, x) G [0,T***] x f2, 

where (f> r (x) = 4>(f) is a cut-off function satisfying <^> G C^°(Bi) and = 1 in Bi. Then from 
(13.51). we derive 



eg {p(0r)t + pv ■ V6 r ) - div^Vfl)^ + divv(p(p, 6>i) - p(p, 6> 2 ))^ r = egpOv ■ V<p r 
Hence multiplying this by 8 r and integrating over [0, t] x we obtain that 



p\9 r \ 2 (t)dx + / / K\V6 r \ 2 dxds 



— 2 — 2 — 1 

d t e e p9 r + Vegpv8 r - 6 r (p r dwv(p(p, 6>i) - <9 2 )) dxds + I r (i), (3.7) 



where the remainder term I r (t) satisfies 



tt2 



k9 (V(f) r y + e e 6p6 r v ■ V(j) r 



\I r (t)\ = 

< Cj J p\ee\\e\\e r \\v\\Vcj) T \ + KM*\V(f> r \ A dxds. 
Since eg, k, p G L°°(0, T*„; L°°) and G L°°(0, r«*; Dj), it follows that 



\I r (t)\ 



C 



< C p\9 r \ 2 dxds + / (p\6\ 2 \v\ 2 + 6 )dxds 



JQ r 



< C I I p\9 r \ 2 dxds + % 



t r 







1772, 



|| ^ 2 1| I ^t- I 3 + ll^lliel^rl 3 



ds 



< C J J p\9 r \ 2 dxds + C 

where the positive constant C is independent of r. Therefore, substituting this estimate into 
(13. 7h . using the similar argument in (13. 6p and applying Gronwall's inequality, we conclude that 



0<t<T* 



sup \\ P m)\\Unr) + I HV0(i)||! 2(nr) 



,dt < c. 



Thus, we have pa/9 g L°°(0, T***; L 2 (ft)) and V6> G L 2 (0, T*„; L 2 (ft)). Now, we can estimate 
I r (t) again, 



i r (t)\ < c [\\ve\\ L 2 + \\p\\ 2 L oo\\p2 0\\ L 2\\v\\ L ™ ) \\ve\\ L 2 {n \ Br) ds 



< C 



|V0 



fll|2 



L2(H\i? 5 ) 



ds 



0, as r — > oo. 
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Hence, letting r — > oo in (|3.7p . we could obtain the estimate (|3.6p and = in (0, T***) x £2. 
Similarly, we have u = 0. This completes the proof of the uniqueness. 

Finally, we prove the time-continuity of the solution (p, 9, u) with the regularities (13,20 - 
(HLD. From Lemma EH we have p G C([0, T***]; L 6 fl D 1 fl D 1,9 ). From (JSHJ) and v G 
C([0, T***]; Dq n -D 2 ), we can obtain p 4 G C([0, T***]; L 2 n L 9 ). To show the time-continuity of 
(0,u), we first observe that (6,u) G C([0, T«*]; Z\} n D 1 ' 9 ') n C([0, T„*]; D 2 - weafc), V </ G (3,6). 
From the equations (j2.2H - fl2.3j) . we also observe that ((egp9t)t, (pu t )t) G L 2 (0, TU*; H~ x ). Since 
(egpdt, put) G L 2 (0, T***; iZg), it follows immediately that (egp9 t , put) G C(0, TU*; L 2 ). Hence 
we conclude that for each t G [0, T***], (0,it) = (9(t),u(t)) £ Djfl D 2 is a solution to the el- 
liptic system — div(KVe) = F and L(a,r],u) = G in f2 for some (F,G) G C([0, T***]; L 2 ). Since 
(k, /x, A)(a, 77) G C([0, T***]; CnD'fl D 1 ' 9 ), using the elliptic regularity estimates (|2T4nj) - (|2T43j ). 
we easily show that (0,ii) G C([0, T***]; D 2 ). From this time-continuity of the solution (p,6,u) 
and the estimates (|2.70p - (j2.71j) . we obtain 

-Po(-)llL6 n z)inD 1 ''J + ll(0(*,O - {-),u(t,-) - u (-))\\ D i nD 2 -► 0, as t -» 0. 
We have completed the proof of Lemma 13.11 □ 

4 Existence result for polytropic fluids 

In this section, we consider the following initial boundary value problem for a viscous polytropic 
fluid: 

Pt + div(H = 0, (4.1) 

pu t + pu-Vu + L(p, 9, u) +Vp = pf, in (0, 00) x (4.2) 

eg(p9 t + pu • V0) - div(AtV0) + pdivu = Q(p, 9, u)+ ph + e p p 2 dwu, (4.3) 

(p, 9, u)\ t =o = (p , O , u ) in O, (4.4) 

(0, «) = (0, 0) on (0, 00) x dCl, (4.5) 

(p, 0, u)(t, x) -> (p°°, 0, 0) as |x| -> 00, (4.6) 

where Q(p, 0, «) = Vu + Vu T | 2 + A(p, #)(divu) 2 , L(p, 0, u) = -div(p(p, 0)(Vu + Vu T )) - 

V(A(p,0)divn) and (eg, e p ,p, p, X, k) = (eg, e p ,p, p, A, n)(p,9). 
Proof of Theorem 11.11 

Our proof will be based on the usual iteration argument and results in Lemma 13.11 Let us 
denote 

c = 2 + p°° + ||po - p°°||L6 nD i nD i, 9 + \\(9 ,u )\\ D i nD 2 + || (gx, g 2 ) \\\-z, 

and we choose the positive constants c\, C2, C3, C4 and T4 as in Section [21 depend only on cq. 
Next, let u° G C([0, 00); Dq D D 2 ) n £ 2 (0, 00; D s ) be the solution to the linear parabolic problem 

uit — Aw = in (0, 00) x O and w(0) = uq in f2. 

Then, chose a small time T5 G (0, T4], such that 

sup ||it°(t)|U + / T5 (ll^WIlL + \\u°(t)\\ 2 D2 , q ) dt < Cl , 

0<t<T 5 JO 7 

SUp ||lt ||£)2 < C2- 

0<t<T 5 
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Hence, choosing T = T5 in Section [21 it follows from Lemma [3. II that there exists a unique strong 
solution (p ,e ,u ) to the linearized problem f)2. 1 j) - (|2.6|) with (a,rj,v) replaced by (po,6o,u°), 
which satisfies the regularity estimates (I2.65|) - (l2.7ip with T*** replaced by T5. Similarly, we 
construct approximate solutions (p k ,9 k ,u k ) as follows: assuming that (p k ~ 1 ,9 k ~ 1 ,u k ~ 1 ) was 
defined for k > 2, let (p k , 9 k , u k ) be the unique solution to the linearized problem (|2.ip - (|2.6p with 
(a,r),v) replaced by (p k ~ 1 , 9 k ~ 1 , u k ^ 1 ). Then it also follows from Lemma 13. II that there exists a 
constant C > 1 such that 



sup (\\p k (t) -p°°||L6 nD i njD i,, + \\p k (t)\\ L 2 nLq 



0<t<T 5 



+ sup \\(e k ,u k )(t)\\ Dl D2 + sup \\{y/?0>,VFvt)\\v 

0<t<T 5 " 0<t<T 5 

+ (\\(0t^)(t)f D i + W k ,u k )(t)\\l 2 , q ) dt<C (4.7) 



\\(6 k (t,-)-e (-), U k (t,-)-u (-))\\ Dh <Ctz, (4.8) 

and 

\\p k (t,-)-P0(-)\\ L e <CtK (4.9) 

for all k > 1. Throughout the proof, we denote by C a generic constant depending only on co, 
but independent of k. 

From now on, we show that the sequence (p k ,9 k ,u k ) converges to a solution to the original 
nonlinear problem (|4.1[) - (j4.6f) in a strong sense. Let us define 

pk+l = p k+l _ yk+l = Qk+l _ ^ _ k+1 = uk+1 _ ^ 

Then, from (|4.ip - ()4.3p . we derive equations for functions (~p k+1 ,0 k+1 ,u k+1 ), 

p k+1 + div(p*+ V) + div(/u fc ) = 0, (4.10) 

pk+ i-k+i +p k+i u k . v -fc+i + L{p k }9 k jlI k + i ) 

= p k+1 (f - u k t - u^ 1 • Vu k ) - p k+1 u k ■ Vu k + V(p* - p k+1 ) 
+div((// - p k - l ){Vu k + (Vu k ) T )) 

+V((X k - X k - 1 )dWu k ), (4.11) 

4 (p k+1 t +l + P k+1 u k ■ Vt +1 ) - div(K k V6 k+1 ) 

= div(( K fc - K k ~ 1 )Ve k ) + Q(p k , 9 k , u k ) - Q(p k - 1 ,9 k ~ 1 ,u k ~ 1 ) 
+p k+1 h + p k dwu k ~ 1 - p k+1 dwu k + e k f 1 p k {9 k + u k ~ l ■ V9 k ) 
-e k e p k+1 (9 k + u k ■ V6 k ) + e k (p k+1 ) 2 divu k - eJ^Go^divu* -1 , (4.12) 

where (e^,e k ,p k ,K k ,p k ,X k ) = (e e , e p ,p, n, p, X)(p k , 6 k ). 

First, we consider the case that p°° > 0. Multiplying (|4.1U|) by p k+1 and integrating over £1, 
we obtain (at least formally) that 



and 



d 
dt 



*p k+1 \ 2 dx 
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< C j (|W s ||p fc+1 | 2 + \Vp k \\u k \\p k+1 \ +/|Vn fe ||7j fe+1 |) dx 

< C (||V U fc || w , 1 , q ||^+ 1 ||| 2 + (HV/IU3 + ||/||L-)||Vn fc || L2 ||p fc+1 || L2 ) . 



Hence, by virtue of Young's inequality, we have 



|l|p fc+1 |l! 2 <^(0llp fc+1 Hi 2 + e ||Vn fc ||i 2 , 



(4.13) 



where A k (t) = C||V 



u 



W i, q + e- 1 C(||Vp fe (t)||| 3 + \\p k {t)\\l°°)- From the estimate flUD, we 
have A k (t) G 1^(0, T 5 ) and f*A k (s)ds <C + C t t for all k > 1 and t G [0,T 5 ]. Here we denote by 
C e a generic positive constant depending only on e -1 and the parameters of C, where e G (0, 1) 
is a small number. 

Next, multiplying (|4.1ip by u k+1 and integrating over Q, recalling that 

{p k+1 )t + div(p fc+ V) = in Q, 

and using the technique in the proof of (|2.56p , we have (at least formally) that 



1 d 
2~dt 



dx + p»J\V^dx 



< c 



\p k+1 \(\f\ + |«t I + |« fc_1 • Vu k \)\u k+1 \ + |/ +1 ||?Z /c ||W s ||u fc+l1 



+\p k+1 - p k \\Vu k+1 \ +{\p k - M fc_1 | + |A fc - A fc - 1 |)|Vu fc ||Vn fe+1 | 



dx. 



Then, using Holder, Sobolev and Young's inequalities, it also follows from (14.71) that 
d 



dt 



|| vV+V +1 f L 2+p°\\ Vu 



Jk+l||2 



li 2 



< B k (t) (||p fc+1 ||i 2 + ||?% + \\y^u k+1 \\h + \\^0 k+1 \\l 2 )+C,\\V?0 k \\h 
+j\\vt\\l 2 +e\\Vu k \\l 2 +C J \p k+1 \\u k \\u k+1 \dx, (4.14) 

where B k {t) G L x (0, T 5 ) such that J * B k (s)ds < C + (7 e t for < t < T 5 and fc > 1. 
Finally, multiplying (|4.12p by # , integrating over Q and recalling that 

(p k+1 )t + div(p k+1 u k ) =0 in 0, 

we obtain (at least formally) that 

'e k p k+1 \6 k+1 \ 2 dx + K° J \V9 k+1 \ 2 dx 



ld_ 

2di 



< C I \ \n k - K*- 1 \\VO k \\vd' +1 \ + (|//| + |A fe |) ( \Vu k \ + \Vu 



k-l 



\Vu k \\0 k+1 \ 



+ |/ +1 ||Vn fc - 1 | 2 (|/ - p k - l \ + \\ k - \ k ~ l \) + |p fc+1 ||/i||^ +1 | 

+ |£ fc+1 | p k d\wu k ~ l -p k+1 d\vu k + e k - 1 p k (0 k + u k ~ l ■ V6 k ) - e k eP k+1 {e^ + u k ■ V6 k ) 
+e k {p k+1 f^u k - e k -\p k fdWu k -l + \{d t e k + u k ■ V e k ) p k+1 t +1 \ } 



dx. 
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From (|4.7p . using Holder, Sobolev and Young's inequalities, we have 



II Lknk+W k+1 \\ 2 -L *- ll\7fl fc+1 ll 2 



< D k 



(*) (\\p k+1 \\h + \\p k \\h + \\^o k+1 \\h)+c.\\V?e k \\l 2 

+ C\\y^uXi+C\\Vu k \\l 2 + ^\\vt\\l 2 +C j \-p k+1 \\e k t \\t +1 \dx (4.15) 



where D k (t) G L x (0,T 5 ) such that f*D k (s)ds < C + C e t and fc > 1. 
Therefore, combining (|4.13p - (|4.15p . defining 



we deduce that 

|(v> fc+1 + llP fc lli 2 ) + e\\ve k+1 \\i 2 +M ||V^ +1 ||i 2 

< E k (t) (> fe+1 + ||p fc ||i 2 ) + C^ k + j\\V0 k \\ 2 L2 + 3eC(||Vn fe ||i 2 + [| W" 1 ^) 

+C J \p k+1 \ (e\6 k \\d k+1 \ + |^ fc ||n fc+1 |) dx, (4.16) 

where E k (t) G L^O, T 5 ) such that J Q * E k (s)ds < C + Ci for < t < T 5 and k>2. 

To estimate the last integral term in ()4.16j) . we assume that f2 is an unbounded domain in 
M 3 , and Claim that there exists a large radius i?o > and a small time Tq G (0, T5] independent 
of k, such that 

\p°° < P k+1 (t,x) < 4p°°, V (t,x) G [0,T 6 ] x (n\B Ro ). (4.17) 
From (|4.17|) . we can estimate the integral term in (|4,16p as follows: for < t < Tq, 



C [ \p k+1 \ fe|^||/ +1 | + l«tll« fe+ ' 



dx 



< C ( W^tWh + W^tWh) W-P k+I \\h + \ ( £ ||V^ +1 ||i 2 + /x°||Vn^||| 2 



and 



C [ \f +1 \Ue k \\e k+1 \ + \ u k \\u k+1 \)dx 



ilTffe+li 



< C\\p k+1 \\ L2 (e\\9*\\ L 40 + \\» + ||^|| L 6||n fc+1 || L a 

< c|l ^ + > 2 f e ||v^|| i2 ||7^^ +1 ||I 2 ||v^ +1 " 2 



(p 1 

iiivu fc+i ni 2 



2 

L 2 



+ ||Vn t fc || L2 ||v / ^ T n fc+1 ||! 2 ||Vu fc+1 " f 



< C(p°°) (J|W t fc ||| 2 + IIV^H 2 :, + l) + J (e||V^ +1 ||i 2 +/ ,°||Vn fc+1 ||i 2 

where the constant C(p°°) depend also on p°°. Therefore, substituting these estimates into 
([4.16p . we obtain 

j t (^ k+l + \\p k \\h) + |l|V^ +1 ||i 2 + ^llVn^Hl, 
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< F k (t) (V fc+1 + \\p k \\ 2 L 2) + C^ k + j||Vr||| 2 + 3eC(\\Vu k \\ 2 L2 + || VU^" 1 1|| 2 ) (4.18) 

where F k G L^O, T 6 ), such that J* F k (s)ds < C(p°°) + C e t for fc > 2 and < i < T 6 . It is easy 
to see that this estimate holds also for a bounded domain f2 since we can choose a sufficiently 
large Ro satisfying Q C Br . 

Now, we can prove (liTfj) as follows. Since p - p°° G L 6 n D 1 n D 1 ' 9 C , where the 
space Co consists of continuous functions on f2 vanishing at infinity, we can choose a large radius 
Ro > 1 independent of fc, so that 

< po(x) < 2p°°, Viefl\% (4.19) 

Z 2 

From Lemma |2.H we have 

p fc+1 (t,x) = p (U k+1 (0,t,x))exp 1^- £ divu k (s,U k+1 (s,t,x))ds^ , (4.20) 

where C/ fc+1 (t, s, x) is the solution to the initial value problem 

§- t U k+1 (t,s,x) = u k (t,U k+1 (t,s,x)), t€[0,T 5 ], 
C/ fc+1 (s,s,x) = rr, s G [0,T 5 ], x G J2. 

Moreover, in view of (14. 7h . we deduce that 

/' 

JO 

and 

j[/ fe+1 (0, t, x)-x\ = \U k+1 (0, t, x) - U k+1 (t, t, x)\ 

< [ \u k (T,U k+1 (T,t, X ))\T <Ct, 

Jo 

for all (t,x) G [0, T5] x Q. Letting Tg = mm(T^,C~ 2 (ln2) 2 , ^C~ l Ro) is a small positive time 
dependent only on T5 and the parameters C, we obtain 



\dWu k (s,U k+ \s,t,x))\ds < I ||Vu fc (s)|| L oo(is < Cti 



[ \divu k (s,U k+1 (s,t,x))\ds < In 2, and |C/ fc+1 (0, t, x) - x| < 

JO 



2 ' 



for all (t, x) G [0,Tg] x f2. In particular, it follows that if < t < Tq and x G Q\Br , then 

U k+1 (0,t,x) G Hence the desired result (liTfjl follows immediately from (|4Tjj^ - (j330|) . 

2 

Now, recalling that ^ fc+1 (0) = and using Gronwall's inequality, we deduce from (|4.18j) that 



tf fc+1 + ll?W ! t ~\\ve k+l \\%ds+ [%\\Vu k +X>ds 



ft 

< C e I Tp k (s)ds 




+ ( j f* HV^Iliadfl + 3eC j\\\Vu k \\ 2 L2 + HV^- 1 !!^)^ exp (c(p°°) + C e t 



< C e I ip k (s)ds 
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+Ce (t J* HV^II^dfl + fj j\\\Vu k \\l 2 + IIVB*- 1 !!^)^ exp (c(p°°) + C e t) , 

for k > 2 and < t < Tq. Hence choosing small constants e > satisfying C e exp(C (p 00 )) = |, 
and T7 G (0, Tq] satisfying exp((7 e T7) = 2, we easily deduce that 

^ k+1 + \\p k \\L*+ Al|W fe+1 ||i 2 d S + f ^\\Vu k +X*ds 



< c jj k {s)ds + \{^f liv/||| 2 ^ + ^ j\\\vu k \\l 2 + HVu*- 1 !!!^) . 

for k > 2 and < i < T?. Using similar arguments as in the proof of (I4,13p - (l4.15p . we can obtain 
sup 0<t<T ip 2 (t) < C easily. Fixing a large K > 1, summing this over 2 < k < K and using 
Gronwall's inequality, we obtain 

K K ,. T7 

u k+X \\ 2 L ,\dl 



Y ™p v fc+1 (t) + E/ 7 ( e ii v ^ +1 i^+^°ii v 



Thus, we have 

Y sup i> k +\t) + Y [ * (e\\vt +1 \\l 2 +p \\Vu k+1 \\ 2 L2 )dt<C<oo, 
where T* = T7. 

Therefore, we conclude that the full sequence (p k ,9 k ,u k ) converges to a limit (p,9,u) in the 
following strong sense: 

j pk_ p i^ p _ p i in L-(0,T,;L 2 ) 

\ (0* u*)->(0, u ) in L 2 (0,T,;Di). l4 " Zij 

It is easy to show that the limit (/?, 0, is a weak solution to the original nonlinear problem 
(I4.1l) - (l4.6j) . Furthermore, it follows from (|4.7h that (p,6,u) satisfies the following regularity 
estimate: 

sup KVpOt,vp»H)(t)\\v+ r (m,ut)(t)f D i + \mu)f mv )dt 

0<t<T, JO J 

+ SUp (||p(t)-p°°||L6 nI?W , ? + ||pt(t)|| i2nW + ||(0,«)(t)|| c l nI32 ) <C. 
0<t<T, v u / 

Then adapting same arguments in the Lemma [3. 11 we can easily prove the time-continuity of the 
solution (p, 6, u) and 

\\p{t,-) -Po(-)IU 6 nDini3i.9 + ||(0(*,O - o (-),u(t,-) - u (-))\\ D i nD 2 -►(), as t -► 0. 

This proves the existence of a strong solution. Now we prove the uniqueness of the strong 
solutions. Let (pi,0i,iii) and (^2,^2,^2) be two strong solutions to the problem (|4.ip - (|4.6p 
with the regularities (|1.17j) - (|1.19j) and we denote by (p, 9,u) their difference. Then following the 
argument used to derive (|4.18p . we easily deduce that 

j t (\\P\\h + \Weo{pi,Oi)piWi? + \\yfp^\\h) + y Pllo* + £||u|| D i 
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< G(t) (||p||| 2 + \Weg{ Pl ,e X ) Pl e\\ 2 L2 + Hv^tlllia) 



where G(t) G L^O,^). Therefore, using Gronwall's inequality, we conclude that p = 8 = u = 
in (0, T***) x S7. This completes the proof of the existence and uniqueness part for the case that 
p°° > 0. 

Now we consider the case that p°° = 0. To prove the convergence in this case, we need to 
modify slightly previous arguments. First, multiplying (I4.10|) by sgn(7J fc+1 )|p fc+1 | 3 and integrating 
over fi, we obtain (at least formally) that 

\-p k+l \ldx < C J\Vu k \\f +1 \l + (\Vp k \\u k \+p k \Vu k \)\f +1 \Ux 

< C\\Vu k \\ wl , q \\p k+1 \\K +C|| /0 fc || L6nD i||Vn fc || L2 ||p fc+1 ^ 



, 3 ■ 

L2 Z/2 



l 

Hence, multiplying this by ||7J fc+1 || 2 3 , we obtain 



d 



^ll/ +1 H!3 <A k (t)\\p k+1 \\li +e||Vn fc ||i 2 , (4.22) 



where A* (t) = C\\Vu k (t)\\ w i, q +e' 1 C\\p k {t)\\ 2 D ,. The uniform bound gZj) implies that f*A k (s)ds < 
C + C t t for all k > 1 and t G [0, T5]. In a similar manner, we can also show that 

J t \lP k+1 \\h < BHt)\\p k+1 \\h + 4Vu k \\h, (4-23) 



where B k (t) G L x (0, T 5 ) such that J * B k (s)ds < C + C e t for all fc > 1 and t G [0, T 5 ]. 

Next, multiplying (0 
proof of (|4.15p . we have 



Next, multiplying (|4,12p by # and integrating over $7, using similar arguments as in the 



, /1!ue y^ +1 ||l 2 + . ||v^ +1 ||! 2 



< (H/ +1 Hl §nL2 + WtWh + ll^^ +1 Hi 2 ) +C||^n fc ||i 2 

+^11^5*11^ + ^\\vt\\l 2 + C||Vn fe ||i 2 , VA; > 1, (4.24) 



where D k (t) G L 1 (0, T 5 ) such that J^ 5 D k (s)ds <C + C e t. 

Finally, from (I4.11|) . using similar arguments as in the proof of (I4.14|) . we easily deduce that 



< E k (t) (\\p k+1 \\ 2 dnL , + l^lli* + Wy/P^^Wb + W^^Wh) 

+C e \\^t\\l, + ^\\vt\\ 2 L2 +e\\Vu k \\l 2 , (4.25) 

o 

where E k (t) G L l {Q,T h ) such that E k (s)ds < C + C e t for k > 1. 

From (|4.22p - (|4.25j) . using same arguments in the proof of (14. 211) . we can show that there exists 
a small time T* G (0, T5], such that the sequence (p k , 9 k , u k ) also converges to a limit (p, 9, u) in 
the sense of 

p k -p 1 ^p-p 1 in L°°(0,T*;L 2 nLl) 



in L 2 (0,T*;D^). 
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Then adapting previous arguments for the case p°° > 0, we can easily prove that (p, 6,u) is a 
unique solution to the problem (|4.ip - (|4.6p with the regularities (|1.17p - (|1.20p . 

It remains to prove the blow-up criterion (jl.21jl . To prove this, suppose that T* < oo, and 
let us introduce define by 

= l + ||^ t (t,.)|| i 2+ /*||V0t(a,-)|li2dfl 

J o 

+ \\p(t, •) - P°°(0IU«nDW.« + \\0(t, OIId W>6 + IK*> OlluJ, 

*(*) = IIPt^Oll^nw + II^WtCt.OllLa + IK^u)^-)!!^ 

+ r|k(r,-)||^ + ||(0,n)(r,-)||i 2 , 8 dr, 
Jo 

and J(i) = + for < i < T*. Let r be a fixed time in (0,T*). Then (p,0,u) is a strong 
solution to the equations P~1 > ([4U |) in [r,T*) x tt, which satisfy (fTTT7]) - (fTTT9]) . 
If (jl.2ip is not hold. We can find a constant K such that 

J(t) + sup $(t) < if. (4.26) 

0<t<T* 

Hence, following similar arguments as in Section [2] with T = T* , we can prove the analogues of 
(1237)1 . (g3HD and (f23TT) : for each i G (t,T*), 

ikvm)(*,oii! 2 + fwu^-^ds 

Jo 

< C K (1 + T*)+ [ [C K (l + \\Ve t \\ 2 L2 ){\\uf D2 + \\y/put\\h) + \\u\\ 2 D%II ]ds, (4.27) 
•/ o 

IKt,0lb» <C K (l + \\^pu t \\ L 2), (4.28) 

and 

IN*, Ollk, < C K (e)(l + H/lli, + IIVP^II^) + ellVutH^. (4.29) 
Combining (|4T27|) - (|4T29|) . we have 

ll(VM)(*,OII| 2 + /'ik^OHL^ 

Jo 

< C K {l + T*) + C K ! (l + \\V6 t \\ 2 L2 )\\^u t \\ 2 L2 ds,Vte(T,T*), 
Jo 

and using the Gronwall's inequality, we obtain that for each t £ (r, T*), 

\\(Vpn)(t, Oll£ 2 + f IK( Sj -)ll V s < ^(l + n- 

From this estimate and (|4.28p - (|4.29p . we have 

IK*. Ob* + f \\u(s r )\\ 2 D2 , q ds<C K (l + T*), Vte(r,T*). 

JO 

And using the estimates (12321) . (l2T4H - (l2^3l ). we obtain 

l|Pt(*i OllL 2 nL<? < Cjf) 
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\\6(t,-)\\ D 2+ f \\6{s,-)f D2 , q ds<C K {l+T*), 
Jo 



and 

J(t) <C K {l + T*), (4.30) 

for all t G (t, T*). It is a contradiction, because the maximality of T* implies that J{t) — > oo as 
t — > T*. So (|1.21|) must be hold. This completes the proof of Theorem II. 1[ 

Remark 4.1. By similar arguments in Lemma 13.11 using the cut-off function, we could prove 
that p k , p G L°°(0,r t ;L 2 ) or p k , p G L°°(0, T*; L§ n L 2 ), and VP^ fe , \/p^ & > x/Pi^ G 
L°°(0,:z;;L 2 ) rigorously. 

5 Some regularity results on elliptic system 

In this section, we derive some regularity estimates for the so-called Lame system: 

P(p, A, u) = -div(p(Vu + Vu T )) - V(Adivu) = F, (5.1) 
and the elliptic system: 

- div(KVe) = G, (5.2) 
where is a bounded or unbounded domain in M 3 , (k, p, A)(x) satisfy 

p{x) > pP > 0, 3A(x) + 2p(x) > 0, k(x) > k° > 0, Vx G fl (5.3) 

The result in Theorem 11.11 rely crucially on these estimates. 

Lemma 5.1. Assume that SI is a bounded domain in M 3 with smooth boundary, and let u G 
Dq(SI) n D 1,q (Q) be a weak solution to the system i5.1]) . where 1 < q < oo. If F £ L q (ft), 
{p,X) G C0(fi) with (5 G (0,1), and {\V p\ + |VA|)|Vu| G L«(Q), i/ten u G D 2 ^(Q) and 

IMb».«(n) < C (||F|| L9(n) + ||V«|| X9( n) + ||(|V/x| + |VA|)|Vu||| L , (n) ) , (5.4) 
where C = C(q,pP,U, Will, X)\\cP(n)) ^ s a positive constant. 

Proof. First, we consider the domain = O Pi Br(xq) CC O,. We choose a function \o £ 
C^(Qr) such that Xo|o fl = 1 and < xo < 1- Letting = Xou, from (|5.ip . we have 

/i(xo)Au - (A(x ) + //(x ))Vdivu = Fi + F 2 , (5.5) 



where 



Xo 



F + Vp- (Vtt + Vu T ) + VAdivu - /x(2(Vxo • V)u + Axqu 



-(X + p) (Vxodivn + Vu ■ Vxo + u ■ VVxo) (5.6) 



and 



F 2 = (p(x) - p(x ))Av + (A(x) + p(x) - A(x ) - ^(x ))Vdivu . (5.7) 
From a well-known elliptic theory due to S. Agmon, A. Doughs and L. Nirenberg in [I], we have 

INHx^n) < Ci (||Fi|| L9 (n) + ||F 2 ||i9(n)) (5.8) 
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where C\ = Ci(g, Q, A)|| c /5(q)). Since /i,Ae C^, We have 

H^IU^n) < C 2 sup | a; - Zol^lNlb^Cn), 

x£B R 



and 



|vo||r)2,g(n) < Ci ll-^i \\li(q) + C2 sup |x - aco 1^ II «o || z?2,« r«) > ( 5 - 9 ) 



where C2 = C2{\\(n, ty\\c0(n))- We can choose Rq > satisfying C1C2RQ < §, and obtain 
IMIz> 2 >?(n H ) - ll u o||£)2,?(n) < 2Ci ||Fi || i9 

< C3(||F|| L< (n) + ||V«|U, ( n) + ||(|VM| + |VA|)|V«||U 9( n)) (5.10) 

where C3 = Cs(q, fjP, A) Hc^fn)) an d f° r all < i? < i?o- Here, we use the Poincare's 

inequality 

||w|U?(n) < C||Vn|| i? (n), 

since u\qq = 0. Second, we assume the fl\ = £>i(0) C\M^_ and consider the domain £2 r = B r HM^. 
We choose a function xi £ Cfi°(B r ) such that xilfir = 1 an< l < xi < 1- Letting ui = Xi 1 ", 
using similar arguments in the proof of (|5.10p . we have that there exists a constant r$ > 0, and 

IM| D a.«(n 5 ) < C4 (H^IU^n) + l|V«||x 9 (n) + 11(1 V/x| + |VA|)|Vu||| L , (n) ) (5.11) 

where C4 = 64(5, //°, Ox, !!(/•*> A) [Ic^fn)) an d for all < r < ro- 

For the general bounded domain O, since dO, G C°°, from (|5.1ip . we easily obtain that for 
each xq G <9f2, there exists r XQ > such that 

\\u\\D^(Br(x )nn) < C 5 (||-P|U 3( n) + ||V«|U <( n) + ||(|V/i| + |VA|)|V«||| M(n) ) (5.12) 

where C5 = C$(q, fjP, ||(/z, A)|| C 3(n)) and for all < r < r XQ . Since Q is a bounded domain, 
using the Finite Covering Theorem, we have 

Il«llu3.«(n) < C 6 {\\F\\ Lq(n) + ||V«|| £ , (n) + ||(|Vm| + |VA|)|Vu||| £ , ( n)) 
where C 6 = C 6 (q,fi ,Q, X)\\cP{n))- a 

Lemma 5.2. Let fl R = B R or B R n M.% with R > 1, and u G £>q(£2r) n 6e a weafc 

solution to the system k5.1\) . where 1 < q < 00. If F G L q (Q R ), (fj,, A) G C^(Q R ) mi/i /3 G (0, 1), 
and (|V/x| + |VA|)|V«| G L1(Q R ), then u G D 2 ^{VL R ) and 

\Hd^(q r ) < C (H^IU^na) + l|Vn|| Lg(nfl) + ||(|V/x| + |VA|)|Vu| || £fl(0jI )) , (5.13) 
where C = C(g,//°,Oi, A^l^m^) is a positive constant, independent of R. 

1 ^ 1 

Proof. If we define p = > 3, v G W ' 9 (fii) by t>(x) = u(Rx), [i R {x) = Rp /j,(Rx) and 

-— 1 

A_r(x) = Rp X(Rx) for x G Qi, then we have: 

P(PR, \r,v) = rI +1 F(Rx) = F R for xen v 
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Hence, it follows from Lemma I5TT1 and the fact \\(hr, ^RjWcP^!) ^ [!(/•*) ^) II CP that 
IMId^i) ^ °7 (ll^lli»(ni) + l|V«||L»( ni ) + ll(|Vw?| + \VXR\)\Vv\\\ Lqm ) , 

where the positive constant Cj = Cj(q, yu , Oi, \\{fj,, X)\\cP((i R )) 1S independent of R. Converting 
this back into the unsealed variables, we obtain 

\W\\ D ^(n R ) < CT^t^llFlU^n^+^iVnlU^n^+^^IKIV/il + lVADlVnllU,^)^ 
< C 7 (\\F\\ Lq{QR) + HVuH^n*) + ||(|V/i| + |VA|)|Vu||| M(njl) ) . 

□ 

Lemma 5.3. Let Q, R = B R n Q, where Q is an exterior domain in R 3 with smooth boundary. 
Choose a fixed number R\ > 1 such that £l c C B Rl and assume R > 2R\. Assume that u G 
Dq(Qr) H D 1,q (Q R ) be a weak solution to the system &5. where 1 < q < oo. If F £ L q (£l R ), 
(/x, A) G C^{n R ) with (3 G (0, 1), and (|V/x| + |VA|)|Vu| G L q {Q R ), then u G D 2 > q (Q R ) and 

\\u\\ D 2, q{nR) < C (H-FHi^n*) + \\Vu\\ L <,(n R ) + ||(|V/u| + |VA|)|Vu| \\ L i(n R ) 

+ ||(|VM| + |VA|)H|U s(naHi )), (5.14) 

where the positive constant C = C(q, /jP , R\, ||(^, A)!!^^ \) is independent of R. 



Proof. Choosing a cut-off function <f> G C£°(R ) such that <f> = 1 in B Rl and <j> = in -B^fti' we 
define 

v = <j)u and to = (1 — 0)tt := t/>tt. 
First, observing that v G TVq ,9 (^2i?i) ! we deduce from Lemma 15. II and Poincare's inequality 

IMl£8(n 2Hl ) ^ c||v-u|| L9(n2fli) , 

(since n|ao = 0), that 

\Md^(q 2Ri ) 

< C (\\P(fi,\,v)\\ Lq{n2Ri) + \\Vv\\ Lq{n2Ri) + ||(|V/i| + IVADIVwHIm^)) 

< C (ll^llw^) + l|V«|U 9( Q 2Kl) + ||(|V/i| + |VA|)(|V«| + \u\)\\ Lq{n2Ri) ) . (5.15) 

Next, to estimate w, we observe that w = in .Brj and io G -Dp H D 1,q (B R ). Then it follows 
from Lemma 15.21 and Poincare's inequality that 

\\ w \\D2-i{n R ) 

< C(\\P(»,\,w)\\ Lq(nR) + \\Vw\\ Lq(nR) + ||(|V//| + |VA|)|VHIIw(n«)) 

< C {\\F\\ Lq(nR) + ||Vu|U 9(0h ) + ||(|Vm| + |VA|)|Vn||| L9(nii) 

+ ||(|VM| + |VA|)H|| w(najli) ). (5.16) 

Combining (|5.15|) and (|5.16p . we finish the proof. □ 

Since the constant C in Lemma I5.2H5.3I is independent of R, using the domain expansion 
technique as in [15] and the fact that the solution to the system (|5.ip is unique in Dq(£Y), we 
easily prove 
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Lemma 5.4. Let Q be the whole space M 3 , the half space or an exterior domain in M. 3 with 
smooth boundary. Assume that u G Dq(Q) n D 1,q (Q) be a weak solution to the system l\5.1\) . 
where 1< q < oo. If F G L^{Q), (/z,A) G C^(fi) mtt /3 G (0,1), and (|V//| + |VA|)|Vu| G L 9 (0), 
i/zen u G D 2,<? (r2) satisfying 

\\u\\ D ^(n) < C {\\F\\ Lq(n) + ||V«|| i9(n) + ||(|V/x| + |VA|)|V«||| w( n)) , (5.17) 
when Q is the whole space M 3 or the half space M 3 ^ ; and 

\W\\v^(Q) < C {\\F\\ Lg{n) + ||V«|U« ( n) + ll(|V/x| + |VA|)|Vu||| £ , (n) 

+ ||(|V/i| + |VA|)H|| L , (naHl )), (5.18) 

when fi is are exterior domain in M 3 with smooth boundary, and J7 C C -B_r 1 with R% > 1, where 
C = C(q,^,\\^,X)\\ c0{Q) ). 

Lemma 5.5. Zei fee £/ie whole space M 3 , the half space M$, an exterior domain in M 3 or a 
bounded domain in M 3 with smooth boundary. Assume F G L 2 (Q), (fj,, A) G C 3 n D 1,p (f2) with 
p > 3, /3 G (0, 1), and u G Dg(f2) D D 2 (r2) 6e a weafc solution to the system \5.1\) . then we have 

\H\d\q) < C MI^II^(Q) + ||V«|| La( n) + ll|V/x| + |VA||||^ n) |M| D i (n) ) , (5.19) 

where C = C(p,fj,°, ||(M; ^)llc 3 (n)) ^ s positive constant. 

Proof. First, we consider the case that f2 is an exterior domain in M 3 with smooth boundary. 
Using Holder, Galiardo-Nirenberg and Young's inequalities, we obtain 

CIKIV^I + |VA|)|Vn]|U 2(n) 
< C|||V/i| + |VA||| LP(n) ||Vn|| ^, 

p-3 3 



< C|||VH + |VA||| iP(n) ||n||^ (n) |M|| )2{n) 



< C|||VM| + |VA|||^ ( 3 n) ||«|| I? i (n) + >| D2(n)l (5.20) 



1 



and 



||(|V^| + |VA|)|n||| L2(Q2fli) 

< C|||V/x| + |VA||| iP(n )H » 

< C\\\Vfi\ + \V\\\\ LP{n2R Ju\\ L(i{n2Ri) 

< C\\\Vfi\ + \VX\\\ LP{Q) \\u\\ Dh{n) . (5.21) 
From Lemma 15.41 and the estimates (|5,20p - (|5.2ip . we have 

h\\n*(ci) < C (\\F\\ L2{n) + ||Vn|| i2(n) + |||V M | + |VA| ||£» n) \\u\\ Dh{n) ) . (5.22) 

By similar arguments, we can easily obtain the result in the case that £1 is the whole space M 3 , 
the half space or a bounded domain in M 3 with smooth boundary. □ 
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Lemma 5.6. Let Vt be the whole space R , the half space R+ ; an exterior domain in R or a 
bounded domain in R 3 with smooth boundary. Assume F G L P (Q), (/U, A) G C@ n D 1,p (r2) and 
u G Dl(£l)P\D 1 ' p {Q)P\D 2 - p (VL) be a weak solution to the system ( f5.il) . where p > 3 and /3 G (0, 1), 
then we have 



\ u \\D 2 'P(n) 



<C \\F 



\LP 



5p— 6 
2p-6 



(n) + l|Vn|| L p (n) + |||V/x| + iVAIII^IInll^i^ 



(5.23) 



where C = C(p, fj, , \\(/jL, \)\\cP(fi\) * s a positive constant. 

Proof. First, we consider the case that Q. is an exterior domain in R 3 with smooth boundary. 
Using Holder, Galiardo-Nirenberg and Young's inequalities, we obtain 



and 



C 

< c 

< c 

< c 

c 

< c 

< c 
c 



(|Vm| + |VA|)|Vu||| £P( n) 
|V/x| + |VA||| LP( Q)||Vn||ioo(n) 



|V/x| + |VA||| LP 



2p-6 
1 5p— 6 



3p 
I 5p-6 



(f^)ll'"ll jD l(Q)lFllD2 i P(n) 
5p-6 2 

|V/i| + |VA|||^(n)||n|| D i (n) + ^||u||2ja.p(n), 



(|V^| + |VA|)|n||| iP(Q2fli) 

|VmI + |VA||| iP(f7Mi) ||n|| L ^ (C2iii) 



(5.24) 



|Vyu| + |VA||Up 



4p-(i 
1 5p-6 



1 5p-6 



(n)iFll£,i( n )|p|i£,2,p( n ) 



|VH + |VA||| 



4p— 6 



+ 



1 



lLf(n)ll U llBi(C) + ^ll«ll2Ja.P(n)- 
From Lemma 15.41 and the estimates (|5.24p - (|5.25p . we have 



D2(n) < C \\F\\ L 2 {n) + ||Vu|| L2(n) + |||V/i| + |VA| 



5p-6 
I 2p-6 



n hi 



DJ(0) 



(5.25) 



(5.26) 



By similar arguments, we can easily obtain the result in the case that f2 is the whole space R 3 , 
the half space R^_ or a bounded domain in R 3 with smooth boundary. □ 

Using similar arguments in Lemmas 15. 5115.61 we can easily obtain: 

Lemma 5.7. Let Q, be the whole space R 3 ; the half space Rj_, an exterior domain in R 3 or a 
bounded domain in R 3 with smooth boundary. Assume G G L 2 (Q), k G C@ n -D 1,p (f2) with p > 3, 
j3 G (0, 1), and e G Dq(Q) D D 2 (f2) 6e a weafc solution to the system \5. then we have 

(5.27) 



ll e llD2(n) < C (11^11x2(0) + ||Ve|| L 2 (Q ) + ||V>t||£ p(n) ||e|| D i (n) 
where C = C(p, k°, ||k||c?/3(q)) wi a positive constant. 

Lemma 5.8. Let Q, be the whole space R 3 , the half space R 3 . , an exterior domain in R 3 u>ii/i 
smooth boundary or a bounded domain in R 3 with smooth boundary. Assume G G L P (S7) ; k G 
qP n L> 1 ' p (^) and e G D^(Q) n D^O) n D 2 - p (VL) be a weak solution to the system i TO) wit/i 
p > 3 and /3 G (0, 1), i/ien we /iaue 

5p-6 



IMId 2 ^^) < C [\\G 
where C = C(p,K°, \\K\\cP(n)) * s a positive constant 



Lp(U) + ||Ve|| £P(n) + HVKll^^Hell^i^) 



(5.28) 
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Remark 5.1. We can remove the assumptions in Lemma [5 . 5H5 . 81 that (u,e) G D 2 or D 2,q . Indeed, 
we can use the usual iteration argument to prove that (u, e) G D 2 or D 2 ' q rigorously. 



6 Blow-up 

In this section, we show the blow-up of the smooth solution to the compressible Navier-Stokes 
equations in W n (n > 1) when the initial density has compactly support and the initial total 
momentum is nonzero. We denote by X(a,t) the particle trajectory starting at a when t = 0, 
that is 

—X(a, t) = u(X(a, t),t) and X(a, 0) = a. 

We set £1(0) = supppo and £l(t) = {x = X(a,t)\a G £1(0)}. From the transport equation (|1.24|) . 
one can easily show that suppp(x, t) = £l(t) and hence from the equation of state (|1.27p that 

p(x, t) = 6(x, t) = e(x, t) = 0, if x G £l(t) c . 

Therefore, from the equation (|1.26|) . we observe that 

^^-\Vu + Vn T | 2 + A(0, 0)(divn) 2 = 0, if x G £l(tf . 

The following lemma was shown in [H [22], when fi, A are constants. Therefore, the proof is 
very similar to that in [H [22], and is just modified to our case that (/x, A) satisfies (jl.28|) . so we 
omit the detail. 

Lemma 6.1. We assume that (fi, A) satisfies HTM) , and (p,S,u) G C 1 ([0, T];H k ), k > [|] + 2, 
is a solution to the Cauchy problem fil.24\) - [l-26}) . Then 

u(x, t) = in x G £l(t) c . 

Moreover, £l(t) = £1(0) for allO<t<T. 

Proof of Theorem 11.31 

We assume that £l(t) = suppp(-,t) is contained in a ball B r ^ t y Multiplying x% to (I1.24D and 
integrating it over M. n , we get the identity 



d 
dt 



pxidx = / puidx. (6.1) 



If we integrate f)1.25|) over M. n , then we also obtain the identity 



^ / pudx = 0. (6.2) 



dt 

Integrating (|6.1[) and ()6.2p over [0,t], respectively, we obtain the following identities 

pxidx = / poXidx + / / puidxds, (6-3) 

JR™ Jo jR n 



pudx = / poUodx. (6-4) 
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From (ESI)- (ED, we get 



/ pxidx = I p Xidx + t pou 0i dx = mu + m 2 ji. (6.5) 

JR n JR n JR n 

On the other hand, since Q(t) C B r M, we can estimate the bound of the L.H.S. of (|6.5p as 
follows, 



/ pxidx 




/ pxidx 




/ pxidx 






Jn(t) 




J\x\<r(t) 



< r(t) / pdx = r(t) / podx = r(t)mQ. (6.6) 

From Lemma 16. 1\ we have Q(t) = O(0) for t £ [0, T] and hence r(t) can be chosen to be r^. 
Therefore, we have 

m r > mu + m 2 it > -m r , i = 1, . . . , n. 
This completes the proof of the theorem 11.31 
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